Accepted for publication in the Astronomical Journal. 
Preprint typeset using I4TgX style emulateapj v. 08/22/09 



ABSORPTION-LINE PROBES OF THE PREVALENCE AND PROPERTIES OF OUTFLOWS 
IN PRESENT-DAY STAR-FORMING GALAXIES 

Yan-Mei Chen', Christy A. Tremonti^, Timothy M. HECKMAN^ Guinevere Kauffmann', Benjamin J. Weiner'', Jarle 

Brinchmann', Jing Wang'-* 

Accepted for publication in the Astronomical Journal. 



o 

^3 

\o 

(N 
< 

6 

(N 
> 

in 

rn 
O 

o 



ABSTRACT 

We analyze star forming galaxies drawn from SDSS DR7 to show how the interstellar medium (ISM) Na 
I AA5890, 5896 (Na D) absorption lines depend on galaxy physical properties, and to look for evidence of 
galactic winds. We combine the spectra of galaxies with similar geometry /physical parameters to create com- 
posite spectra with signal-to-noise ^ 300. The stellar continuum is modeled using stellar population synthesis 
models, and the continuum-normalized spectrum is fit with two Na I absorption components. We find that: (1) 
ISM Na D absorption lines with equivalent widths EW > 0.8 A are only prevalent in disk galaxies with specific 
properties - large extinction (Ay), high star formation rates (SFR), high star formation rate per unit area (Ssfr), 
or high stellar mass (M*). (2) the ISM Na D absorption lines can be separated into two components: a quiescent 
disk-like component at the galaxy systemic velocity and an outflow component; (3) the disk-like component is 
much stronger in the edge-on systems, and the outflow component covers a wide angle but is stronger within 
60° of the disk rotation axis; (4) the EW and covering factor of the disk component correlate strongly with 
dust attenuation, highlighting the importance that dust shielding may play in the survival of Na I. (5) The EW 
of the outflow component depends primarily on Ssfr and secondarily on Ay; (6) the outflow velocity varies 
from ^ 120 to 160 km s"' but shows little hint of a correlation with galaxy physical properties over the modest 
dynamic range that our sample probes (1.2 dex in log Ssfr and 1 dex in logM*). 

Subject headings: galaxies : evolution — galaxies : star formation 



1. INTRODUCTION 

Galactic-scale gaseous outflows ('galactic winds') are 
known to be ubiquitous i n very actively star forming galax- 
ies at all cosmic epochs (Heckman et al. 119901 ; iPettini et al 



1^00; Shapleyetal. 2003; Menard et arf lfoogT IWemer et al 
12009). Galactic winds play an vital role in the evolution of 
galaxies and the intergalactic medi um (I GM). The "baryon 
deficit" in the Galaxy (e.g., Silkj l2003h indicates gas re- 
moval by outflows during past active episodes. The mass- 
metallicity and effective yield relations observed in local 
galaxies suggest that galactic winds transport highl y metal- 
enriched gas out of galaxies and into flie IGM (e.g.. lGarnettl 
l2002t ITremontietalJi200lt iDalcantonI 12007' ). The ACDM 
model over-predicts the galaxy luminosity function at both 
the low and high luminosity ends if a constant mass-to-light 
ratio is assumed. The most natural way to reconcile the obser- 
vations with theory is to i nvoke feedback processes, including 
supernova, stellar winds (Cole et al.'2000';'Benson et al.l2003l ; 
IStringer et al. 2009; OcDenheimer et al. 2009) and AGN ac- 
tivity (Silk & Rees 1998; Hopkins et al. 2006). However, it 
remains unclear which kind of feedback processes dominate 
as a function of luminosity and comic epoch. 
In actively star forming galaxies, galactic winds are 
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driven by the mechanical energy and momentum imparted 
by stellar winds and supernovae (Chevalier & Clegg 1985|; 
iHeckman et alJ Il990h . Young star clusters create over- 
pressured bubbles of coronal phase gas which expand and 
sweep up shells of ambient ISM until they 'blow-out' of the 
disk into the halo. The collective action of multiple superbub- 
bles drives a weakly collimated bi-polar outflow consisting of 
hot gas and cool entrained clouds. Radiation pressure is also 
likely to play a role in accelerating cool dusty material (e.g., 
Murra v et al.ll2()()5h . In this paper we use the term 'galac- 
tic wind' to describe such outflows without regard for their 
eventual fate — the gas may escape the halo potential well, 
or be recycled back into the disk in a process sometimes re- 
ferred to as 'galactic fountai n' activity (Shapiro & Field 197^ 
lBregmanl[T980l: Ik^[T981 . 

In the last two decades, there have been many attempts 
to directly observe galactic winds in galaxies at z = - 3 
dVeiUeux et alj|2005l) . In the local universe, outflows can be 
detecte d via X-rays which tr ace hot gas (e.g ., Dahle m et alJ 
[1991 iMartin et alJ l200l IStrickland etall i2004), opti- 
cal nebular emission lin es produced by warm gas (e.g., 
iLehnert & Heckmanll996|) and ISM absorption lines (e.g., Na 
I. K I) from cold gas (e.g.. lHeckman 2002) . Blue-shifted Na D 
absorption from the entrained cool gas is frequent ly detected 
in IR-bright starburst galaxies (Heckman et al."2000), LIRGs 
and U LIRGs (Rupke et al. 2002, 2005a b; Martin 2005, 200^ 
iMartin & Bouche 2009), and it is sometimes evident in dwarf 
starbursts (Schwartz & Martin 2004). The velocity of this gas 
correlates weakly with star formation rate (SFR) and galac- 
tic rotation speed, with a factor of ~ 30 change in veloc- 
ity observed over a range of 4 orders of magnitude in SFR 
(v„i„d oc SFR" \ Rup ke et alJl20()5btlMar tin'2005). 

At intermediate redshift, the Mg II AA2796, 2803 ISM 
absorption line shifts into the observed-frame optical. The 
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relatively wide separation o f this doublet makes it a good 
choice for outflow studies. iTremonti et alj ( l2007h detected 
500 - 2000 km s~' outflow velocities in a small sample of 
z ~ 0.5 post-starburst galaxies, and suggested that these out- 
flows cold be the relics of AGN-driven winds. IWeiner et'al] 
(I2OO9.) employed the spectral stacking technique to probe 
Z ~ 1.4 star forming galaxies drawn from the DEEP2 sur- 
vey. This study, which is based on ~ 1400 galaxy spectra, 
demonstrated that blue-shifted Mg II absorption is ubiquitous 
in actively star forming galaxies. 

At high redshift, the observed-frame optical samples the 
rest-frame far ultraviolet, which is rich in strong ISM res- 
onance absorption transitions, but lacking in stellar spectral 
features. Because of the difficulty of measuring the relative 
velocity of the stars and gas, outflows have only been de- 
tected in very luminous Lyman break galaxies (Pettini et al] 
|2()00, 2002) or in composite spectra (IShaplev et al.H2003 ). 

While much of our knowledge about galactic winds comes 
from studies of local star forming galaxies, the work to date 
has been based on relatively small samples of extreme objects 
(dwarf starbursts, ULIRGs). The properties of outflows in 
local normal star forming galaxies are still largely unknown. 
In this paper, we investigate Na D absorption in a sample of 
^150,000 star forming galaxies drawn from the Sloan Digi- 
tal Sky Survey. By stacking the spectra of galaxies selected 
to have similar physical attributes, we obtain very high S/N 
composite spectra. After carefully modeling and dividing out 
the stellar continuum, we are able to probe very low Na D 
Equivalent Widths (EWs) and extend our analysis over a wide 
range in galaxy physical parameters. 

This paper is arranged as follows. In §2, we introduce the 
sample selection criteria used in our study. Our method of 
creating composite spectra and measuring the Na D lines is 
developed in §3. We apply the method to the SDSS sample in 
§4 and summarize our results in §5. We use the cosmological 
parameters //q = 70 km s"' Mpc"', 51m = 0.3 and JIa = 0.7 
throughout this paper 

2. SAMPLE SELECTION 
2.1. The Data 

The Sloan Digital Sky Survey (SDSS; lYork et al.l l2000t) 
spectroscopic galaxy catalogue con tains ^930,000 spectr a 
in its seventh data release (DR7; lAbazajian et al.l |2009|). 
We analyze qbiects drawn from the "Main" galaxy sample 
(IStrauss et al.ll2002h . which have Petrosian r magnitude in the 
range 14.5 <r < 17.7 after correction for foreground galactic 
extinction using the reddening maps of Schlegel et al.l (Il998i) . 
The spectra are taken through 3" diameter fibers and cover a 
wavelength range from 3800 to 9100A with a resolution of 
R ^ 2000 and a dispersion of 69 km s"' pixel"'. Most of the 
galaxies have redshifts in the range z = 0-0.3. We adopt red- 
shifts from the specBS pipeline (lAdelman-McCarthy et al.l 
I2008h . The median redshift error is ~ lO"'' 

As described in ITremonti et al.l (120041) and 
iBrinchmann et al.l (120041) . a stellar population model is 
fit to the continuum of each galaxy spectrum after masking 
out the strong nebular emission lines. The basic assumption 
is that any galaxy star formation history can be approximated 
as a sum of discrete bursts. The library of template spectra is 
composed of single stellar population (SSP) models generated 
using a preliminary version of the population synthesis code 
of Chariot & Bruzual (2010, in prep.) which incorporates 
the MILES empirical stellar library (.Sanchez-Blazquez et all 



I2OO6I) . We refer to these models hereafter as the 'CB08 mod- 
els'. The spectral resolution (2.3 A FWHM), spectral-type 
coverage, flux-calibration accuracy and number of stars in 
the MILES library represent a substantial improvement over 
previous libraries used in population synthesis models. We 
generate template spectra with ten different ages (0.005, 
0.025, 0.1, 0.2, 0.6, 0.9, 1.4, 2.5, 5, 10 Gyrs) and four 
metaflicities (0.004, 0.008, 0.017, 0.04). For each metalHcity, 
the ten template spectra are convolved to the measured 
stellar velocity dispersion of each individual SDSS galaxy, 
and the best fitting model spectrum is constructed from a 
non-negative linear combination of the template spectra, with 
dust attenuation modeled as an additional free parameter 
The metallicity which yields the minimum is selected as 
the final best fit. The results of this fitting procedure and 
measurements of an number of line indices (e.g., Lick_Mgb, 
D4000) are made available in the MPA/JHU catalogue^. The 
catalog also includes the equivalent width (EW) and flux of 
an number of emission and ISM absorption lines (e.g.. Ha, 
H/3, Na D). 

The SDSS im aging data consist o f CCD imaging i n 
u,g,r,i,z bands (Fukugita et al. 1996; Smith et al. 2002), 
taken with a drift scan camera (Gunn et al. 1998) mounted 
on a wid e-field 2. 5-m telesc ope. The SDSS photometric 
pipeline (' Lupton et al.l 120011) fi ts each galaxy image with 
a two-dimensional model of a Ide VaucouleursI |l948) sur- 
face profile and an exponential profile, each convolved with 
the PSF of the image. A seeing-independent axial ratio 
(b/a) is derived by this fitting procedure. The pipeline 
also computes the best linear combination of the exponen- 
tial and de Vaucoule urs models and stores it in a param- 
eter called fracDeV ( Abazaiia n et al.l 120041) . Here we fol- 
low l^dilla & Strauss (2008) in using fracDeV to distin- 
guish disk galaxies (fracDeV < 0.8) from early-type galax- 
ies (fracDeV > 0.8). The axis ratios from the exponential 
and de Vaucouleurs models are consistent with each other 
independent of the morphology of the galaxy. In this work, 
we are only concerned with disk galaxies, and therefore we 
adopt the axis ratios from the exponential fit. We com- 
pute the galaxy inclination, /, from the measured axial ra- 
tio, b/a, and the r-band abso lute magnitude Mi using Table 
8 in iPadilla & StraussI (120081) . This empirical relation takes 
into account the increase in bulge-to-disk ratio with absolute 
magnitude. Typical inclination errors are 10° at / > 40°. At 
/ < 40° galaxy asymmetries can result in over-estimates of 
the inclination. The photometric properties used in this pa- 
per (fracDeV, b/a, u,g, r, i,z magnitudes) are also available in 
catalogs on the MPA/JHU webpage. 

The derived galaxy parameters required in this work in- 
clude stellar mass (M* for the whole galaxies and Mjjber for 
stellar mass within fiber), stellar mass surface density (S^. ), 
V-band dust attenuation (Ay), SFR, SFR surface density 
(Ssfr), and the SFR per unit stellar mass or 'specific SFR' 
(SSFR = SFR/Mfiber)- The stellar mass-to-light ratios are ob- 
tained by comparing urgiz colors of galaxies to a large grid of 
CB08 model colors following the methodology described in 
Salimet al. (2007) to avoid aperture correction issues. This 
approach differs from that of Kauf fmannI (j2003a) who used 
the spectral indices D„(4000) and HSa to constrain the mass- 
to-light ratios. However, the differences between the mass 
estimates of the two methods are very small, typically 0.01 

^ T he MPA/JHU catalog can be downloaded from 
|http://www.mpa -garching. mpg.de/SDSS/DR7 
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dex. 

We calculate Ay, the optical dust attenuation in magnitudes, 
using the Ha/HP ratio. We assume an intrinsic Ha/H/3 ratio 
of 2.87 (Osterbrock & Ferland 2006), a Milky Way attenua- 
tion curve (iCardelli et al.lll989 ). and R = 3.1. The nebulai- 
attuenation is likely to overestimate the attenuation of the stel- 
lar continuum (c.f.. lCalzetti et alj|1994l) . 

We derive the SFR from the dust extinc t ion cor rected Ha 
luminosity using the formula of iKennicuttI (Il998al) . dividing 
by a factor of 1.5 to convert from a Salpeter to Kroupa initial 
mass function (Kroupa 2001). Ssfr is defined as SFR/ttR^ 
and Sm, = ^fiber/i'^^^ where R corresponds to the SDSS 1.5" 
aperture radius corrected for projection effects. 

Except stellar mass, all of these derived parameters are cal- 
culated from quantities measured within the SDSS fiber aper- 
ture and we do not apply aperture corrections. We believe the 
locally determined physical parameters will have the greatest 
influence on the local ISM kinematics traced by Na D; star 
formation occurring in the outer disk of the galaxy is unlikely 
to drive a circum-nuclear outflow. Another exception to our 
policy of using locally derived physical properties is FigurefTTl 
where we use global galaxy SFRs in order to compare with 
results from the literature. The global SFRs are estimated by 
adding the SFR in ferred from the emission lines measured in 
the fiber aperture jBrinchmann et al.ll200 4') and the SFR esti- 
mated from the colors outside the fiber from fits to the 5-band 
photometry. The total SFRs are comparable the global SFRs 
in Brinchmann et al. ( 2004), but for a small subs et of galax- 
ies an aperture correction problem identified by ISalim et al] 
(I2007h has been rectified. 

2.2. Sample Selection 

The criteria used to select the parent sample used in our 
analysis are the following: 

1. redshift of 0.05 < z < 0.18. With 3" diameter fibers, 
the SDSS spectra probe the central 3-12 kpc in this red- 
shift range. At the median redshift of z = 0.09, the fiber 
probes the central 6 kpc. 

2. r-band fracDeV < 0.8. This ensures that we are select- 
ing disk galaxies, which makes the calculation of the 
inclination angle from b/a possible. 

3. D„(4000)< 1.5. This index measures flie 4000 A break, 
and is a good indicator of the galaxy star formation his- 
tory (.Kauffmann.2003 a) . Our cut selects galaxies with 
young stellar po pulations. It is just b elow the thresh- 
old identified bv ' Kauffmaniil (1200 3^ D„(4000)=1.55) 
where a sharp change in galaxy structural properties oc- 
curs. 

4. log([OIII]/H/3) < 0.61/{log([NII]/Ha)- 0.05}+ 1.3. 
We exclude AGNs from our s ample using the ern ission 
Une ratio diagnostics given in lKauffmannl(l2003ch . 

We refer to this sample hereafter as 'Sample A'. It contains 
140,625 galaxies. In §4 we will justify the application of some 
additional cuts to help isolate galaxies with measurable inter- 
stellar Na D. In later sections of the paper, we will refer to 
these sub-samples as Samples B and C. 

In Figure [U we illustrate the distribution of galaxy physical 
properties in Sample A and show how these parameters corre- 
late with one another. As we will discuss in ^4.21 these corre- 



lations make it difficult to distinguish the physical parameters 
that most strongly influence the properties of galactic winds. 

2.3. Interstellar Na D in SDSS Spectra 

Studying gas kinematics using optical absorption line spec- 
troscopy is challenging because of the lack of strong ISM ab- 
sorption lines. An additional complication when using the 
galaxy continuum as the background source, is the the promi- 
nence of stellar absorption features in the continuum. The 
Na I "D" AA5890, 5896 doublet is favored because it is moder- 
ately sttongJsom^tim£s_saturated) in starburst galaxy spectra 
(c.f.. lHeckman et an i2000). and because the stellar contami- 
nation is modest in galaxies dominated by very young stars. 
Other transitions which are sometimes used are Ca II H & 
K, which occur in a region of the spectrum rife with stellar 
absorption and nebular emission lines, and K I AA7665,7699 
which is a relatively weak transition. 

The main drawback of using Na D to probe the ISM of nor- 
mal star forming galaxies is that it is a very prominent fea- 
ture in the the spectra of cool stars; stellar contamination can 
therefore be significant. We discuss this issue further in ^3.21 
where we show that ^80% of the observed Na D absorption 
in our galaxy spectra arises in stellar atmospheres. 

In most previous work, the stellar contribution to the Na I 
line was estimated using the Mg I A A5 1 67 , 5 1 7 3 , 51 84 tr iplet, 
which arises only in stellar atmospheres. iMartinI (120051) find 
EW(Na I) = 1/3 EW(Mg I) based on an analysis of steflar 
spectra. Here we take a somewhat more sophisticated ap- 
proach, and model the galaxy continuum light using stellar 
population synthesis models ( ^3. 2b . We use the Mg I lines 
residuals to gauge the success of our continuum fitting proce- 
dures. 

In our subsequent analysis we will make use of composite 
galaxy spectra. However, we begin by examining individual 
spectra in Sample A with S/N pixel"' > 15. Roughly 19% 
(27,076) of the Sample A galaxies make this cut. For each 
galaxy we measure t he Na D EW and the Lick Mg b index 
dWorthev et alJ[T994h of the observed spectrum and the best 
fit continuum model. In Figure |2] we plot the residual ab- 
sorption (data - model) as a function of nebular attenuation. 
We hypothesize that the Na D excess in the data relative to 
the models is due to ISM absorption. Na I has an ionization 
potential less than Hydrogen, so dus t shielding is expec ted to 
play an important role in its survival dMurrav et alj200'^ . The 
strong positive correlation between the Na D residual and dust 
attenuation is evidence that the excess Na D absorption arises 
in the dusty ISM. The lack of any corresponding correlation in 
the Mg b residual indicates that the starlight is well fit by our 
models even in heavily dust attenuated galaxies. (The small 
positive offset of the median is due to the presence of a weak 
[N I] A5198 emission in one of the continuum bands of the 
Mg I index - see Fig.|5]). 

One notable feature of Figure |2] is tendency for galaxies 
with low dust attenuation to have negative Na D residuals - 
more Na D absorption in the stellar population models than in 
the data. There are three possible explanations for this: (1) the 
S/N of the spectra are not high enough for a robust continuum 
fit; (2) the stellar population models are not robust around the 
Na D region, at least for some stellar populations; (3) Na D is 
primarily in emission in gas and dust-poor galaxies. Through 
the spectral stacking technique, the S/N issue can be eas ily 
resolved. We discuss the other possibilities further in ^3.31 

2.4. Galaxies with Strong ISM Na D Absorption 
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Fig . 1 . — Correlations among the physical parameters of our Sample A galaxies. For panels where the x— and y— axis are the same, histograms are shown. Note 
that all quantities except M* are measured within the SDSS fiber aperture. 



In the previous section we found that the Na D residual 
increases with Ay. In this section, we study how the Na D 
residual evolves with the other galaxy parameters that we are 
interested in: SFR, Ssfr, M^,, Sm,, SSFR. We define galaxies 
with a Na D residual > 0.8 A as strong ISM Na D absorption 
galaxies. This definition is somewhat arbitrary, however, our 
conclusions are not very sensitive to the threshold we adopt. 

Deep images of the edge-on starburst galaxy M 82 in the X- 
ray. Ha. and Sfim-bands (,C) h vama et al. .2002 : Mutchler et alj 
l2007t lEneelbracht et al.ll2006l) " indicate that gas and dust are 
being driven out along the galaxy minor axis in a weakly- 
coUimated bipolar outflow (See Figure O. Imaging of the 
Ha and soft X-ray emission shows that this wind geometry 



is common in starburst galaxies (e.g., iLehnert & Heckmanl 
1996; Strickland et al. 2004). Thus, absorption line studies 
of nearly face-on systems (low /) should probe gas in the out- 
flow, while observations of nearly edge-on systems (high /) 
will not; they may, howeve r, be sensitive to gas in the disk. In 
keeping with this picture, iHeckman et al.l (120001) found that 
there is a high probability (~ 70%) of detecting outflowing 
gas in absorption in starburst galaxies with inclinations less 
than 60°. Based on these results, we split our Sample A into 
two sub-samples with / < 60° and / > 60°. 

Figure H] shows how the fraction of galaxies with Na D 
residuals greater than 0.8 A changes with various galaxy pa- 
rameters for our two samples split by inclination. The frac- 
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Fig. 2. — Residual Na D (left) and Mg h (right) absorption versus dust 
attenuation for Sample A galaxies with S/N > 15. Red dots indicate the 
median. The residual absorption is measured after subtraction of the best fit 
stellar continuum model. Dust attenuation is measured from the Ha and H/3 
nebular lines. The strong positive correlation in the left panel is evidence that 
the excess Na D absorption arises in the dusty ISM. The lack of correlation 
for Mg h is expected since this is an excited photospheric feature that does 
not have a counterpart in the ISM. 




Fig. 3. — The galactic wind in the prototypical starburst galaxy M82. The 
color-composite image was constructed by the Hubble Heritage team from 
B,V,I, and Ha images obtained with the Hubble Space Telescope Advanced 
Camera for Surveys iMutchler et al. 2007). The galaxy has an inchnation of 
i Si 80° . Outflowing gas is visible as a bi-polar cone of Ha emission (red) 
extending along the galaxy minor axis. If this wind geometry is common in 
disk galaxies, then we would expect to observe blueshifted Na D when our 
sightline probes the wind bi-cone at i < 60° . The dusty filaments evident in 
the disk may give rise to Na D absorption at the galaxy systemic velocity. 



tion increases rapidly with all of the galaxy properties, except 
SSFR, for both low (black) and high (red) inclination samples. 
The h ighly inclined galaxies have stronger Na D residuals; in 
^4.11 we will argue that this strong Na D absorption arises in 
the disk. Galaxies with strong ISM Na D absorption make up 
only a small part of the full sample; however, they dominate 
the populations with strong star formation, high stellar mass, 
and high dust attenuation. The blue vertical line in each panel 
marks the place where the fraction of galaxies with / < 60 ° 
and strong Na D reaches 5% and begins to rise steeply. 

3. DATA ANALYSIS 

In this section, we describe how we measure t he pr operties 
of the ISM Na D absorption. As discussed in ^2.31 the low 
S/N of the individual spectra will lead to large uncertainties 
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Fig. 4. — The fraction of strong Na D residual galaxies as a function of dif- 
ferent galaxy parameters. The fraction is measured as F = N(EW > O.S)/^!,;^, 
where ^(EW > 0.8) is the number of galaxies with Na D residual > 0.8 A 
in a given bin and N\yi^ is the total number of galaxies in the bin. The black 
dots represent galaxies in Sample A with < 60° and red dots represent the 
galaxies with (' > 60 ° . The blue vertical line in each panel marks the place 
where the fraction of ! < 60° galaxies with strong Na D residual reaches 5% 
and begins to rise steeply. 

in the measured Na D EW. The median Na D EW error for 
Sample A is ~'1.5 A, which is greater than the strength of the 
interstellar absorption in most galaxies. By stacking (averag- 
ing) galaxy spectra we avoid this S/N issue and are able to use 
samples that are relatively complete. The steps of our analysis 
are as follows: 

• Stack spectra in different bins of galaxy physical pa- 
rameters to create high S/N composite spectra. The pa- 
rameters used in our stacking analysis are galaxy incli- 
nation (0, M*, Sm., SFR, Esfr, SSFR, and Ay. 

• Fit the CB08 stellar population models to the stacked 
spectra generated in the first step. The aim is to separate 
stellar and ISM Na D absorption. Divide each spectrum 
by the best fitting continuum model. 

• Fit the interstellar Na D line in each continuum- 
normalized spectrum. Return the EW, line centroid, 
line width, optical depth, and covering factor 

• Look for correlations between the properties of the gas 
and the galaxy physical properties. 

Each of these steps is described in more detail below. 
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3.1. Stacking galaxy spectra 

The galaxy spectra are first correcte d for foreground Galac - 
tic attenuation using the dust maps of Schlegel etaP OMi), 
transformed from vacuum wavelengths to air, and shifted 
to the restframe using the redshift determined by the SDSS 
specBS pipeline. We then normalize each spectrum to its 
median flux between 5450 to 5550A, where the spectrum is 
free of strong absorption and emission lines. The normal- 
ized restframe spectra in a given parameter bin are averaged 
("stacked") using the following weighting scheme: 



,/comp('^) ■ 



^"^j mask(i, A) x /(/, A) 
X^/Li >nask(i, A) 



(1) 



where the sum is over / = 1,2, ...,n galaxies with normalized 
flux /(/, A). We set mask(i, A) = for bad pixels (identified 
in the SDSS mask array), and mask(i, A) = 1, otherwise. This 
method gives equal weight to all spectra, but excludes bad 
pixels from the composite. 

We adopt an adaptive binning approach to produce com- 
posite spectra in different parameter bins with equal S/N. The 
galaxies are sorted by the physical parameter of interest and 
added to the stack one at time, from lowest to highest. After 
each addition, the S/N of the continuum near Na D is com- 
puted. The process is repeated until a S/N of 300 is reached. 
The value of the physical parameter of interest in each bin is 
computed from the median of the individual measurements. 

3.2. Fitting the Continuum 

The aim of this study is to use the ISM Na D absorption 
to probe galactic winds in star forming galaxies. However, as 
discussed in ^2.31 we can not disregard the fact that interstel- 
lar Na D absorption is superimposed on stellar N a D absorp- 
tion. In previous studies of LlRGs and ULlRGs ( Martinl2005l 
|2006; Rupke et al. 2005a b), stellar Na D absorption was es- 
timated from Mg 1, and deemed to contribute only a small 
fraction (< 10%) to the total line EW. Howev er, in the inves- 
tigation of 32 far-lR-bright starburst galaxies. iHeckman et al.l 
(l2000i) found that > 40% of the galaxies (mostly edge-on sys- 
tems) had strong stellar contamination. It is well-known that 
Na D is strong in the spectra of cool stars, with a peak strength 
in the range from K3 to MO ( Jac oby et al . 1984). These stel- 
lar types are expected to make a significant contribution to the 
current sample since the SDSS spectra are obtained through a 
3" circular fiber aperture that samples the central 3-12 kpc of 
our current sample. In these regions, bulge K-giants are likely 
to make a strong contribution to th e integrated the stellar light 
fijeckman 1980: iBica et al.llT99l]) . 

To separate the stellar and ISM contribution to the Na D 
absorption we model the stellar continuum of each stacked 
spectrum using the CB08 stellar population synthesis models. 
In the Appendix, we show that our results are insensitive to 
the choice of stellar population models. Our continuum fitting 
procedure is nearly the same as that used in fitting individual 
spectra (see ^2. lb . The main difference between the individ- 
ual and stacked spectra is that the stacked spectra have much 
higher S/N. This makes it possible for us to constrain the stel- 
lar contribution using a larger number of SSP templates. We 
therefore fit the stacked spectra with all 40 templates (4 metal- 
licities and 10 ages) simultaneously. This allows the best-fit 
model to consist of a mixture of metallicities, unlike the situ- 
ation for the individual spectra 02.1b . 

The He 1 A5876 nebular emission line on the blue side of 
Na 1 A5890 is broad enough to slightly influence our mea- 



surement of Na D. Since the line presumably originates in 
H II regions behind the Na D absorbing gas, we include it in 
our continuum model. We fit a Gaussian to the blue wing of 
He I emission and add a model of the full line to the stellar 
continuum before normalizing the spectrum. 

Observationally, in local galaxies, galactic winds are a gen- 
eral consequence of high SFR or Ssfr (Lehnert & Heckman 
1996); they are common in galaxies above a Ssfr threshold 
of 0.1 Mq yr"' kpc"' (Heckman 2002). To illustrate our con- 
tinuum and absorption line fitting procedure, we randomly se- 
lect 2000 galaxies with Esfr > O.IM© yr"' kpc"^ from Sam- 
ple A and stack them. Figure |5] shows the resulting stacked 
spectrum (black line) and our stellar continuum fit (red line). 
We highlight the regions near Mg I (left insert drawing) and 
Na I (right insert). The vertical black lines mark the restframe 
wavelength center of the Mg I AA5167, 5173, 5184 triplet, 
[N I] A5198, He I A5876, and the Na I AA5890, 5896 doublet. 
A Gaussian fit to the He I emission line is shown in cyan. 
From this plot, we infer that: (1) the stellar continuum model 
provides a very good fit overall; (2) the stellar contribution of 
the Na D absorption is as high as ^80% on average. Since 
the full Sample A includes galaxies with lower Ssfr and a 
reduced incidence of strong ISM Na D absorption, we expect 
the stellar contribution to the Na D line to be even larger in 
some cases. 

3.3. Measuring Interstellar Na D Absorption Line Profiles 

There are two common methods of measuring absorption 
line profile parameters: (1) fitting simple functional forms 
(such as Gaussian or Voigt profiles) to the intensity as a func- 
tion of wavelength; (2) making more complex models, in 
which the intensity profiles are direct functions of physical pa- 
rameters (velocity, optical depth, and covering fraction; e.g., 
iRupke et al. 2005a). The first solution is mainly based on 
mathematics rather than physics, and it is widely used for 
single lines or unblended doublets or multiplets. The study 
of unblended transitions is much easier than blended ones. 
For the blended doublet or multiplet, we prefer the second 
method since on one hand, the profile shape is readily un- 
derstood in terms of a set of physical parameters, and on the 
other hand, different model assumptions can be tried. For the 
blended Na doublet, we choose to use the second method. We 
briefly summarize Run ke's model below and refer the reader 
to Rup ke et al.l ( I2005al) for more details. 

Previous studies of Na D absorption troughs in ULIRGs 
strongly suggest the assumption of complete continuum cov- 
erage is not met (Martin 2005U200 6; Rupke et alJI200 5ab). 
These studies also found super-thermal Na D fine widths 
(b ^ 300 km s"'), which can arise from multiple fragmented 
shells of cool gas traveling at different velocities (Fuiita et alj 
l2009h . When stacking many hundreds of galaxy spectra, ab- 
sorption lines with different velocities, line widths, optical 
depths and covering factors will be blended. Jenkins. (19861) 
explored this situation extensively using monte-carlo simula- 
tions, and determined that reasonable results can be obtained 
for the ensemble properties of the absorbers provided the dis- 
tribution function of their individual properties is not too un- 
usual. Results from stacked spectra may, in fact, be more ro- 
bust than those from individual galaxies, because the number 
of blended components is larger. 

We adopt the partial covering model for the current study 
and assume a covering factor, C/, which is independent of 
velocity. In this case, the general expression for the intensity 
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Fig . 5 . — An example of a stacked spectrum (black) and its best-fit continuum model (red). The spectrum is a composite of 2000 randomly selected galaxies in 
Sample A with Ssfr > O.IMq yr"' kpc"". The two insert drawings highlight regions near the Mg I (left) and Na I lines (right). The vertical black lines mark the 
restframe wavelength of the Mg I AA5167, 5173, 5184 triplet, the [N I] A5198, the He I A5876 emission line, and the Na I AA5890, 5896 doublet. A Gaussian fit 
to the He I emission line is shown in cyan. 
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Fig. 6. — An example of our Na D line profile fits. The back line in the top panel shows the spectrum from Figure[5]normalized by its best fit continuum model. 
The left panel illustrates our two-component Na D fit: the red line represents the systemic component; the blue line is the outflow (blue-shifted) component; and 
the green line shows the combined fit. For comparison, we also show a single component fit in the right panel. The residuals of are shown in the bottom panels. 
The two component fit does a better job of modeling the shape of the line profile. 
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of a doublet, /(A), assuming a continuum level of unity, is 

/( A) = 1 - Cf + Cfe""'^^^-'"''-^^ . (2) 

where tb(X) and tr(A) are the optical depth of the blue and red 
components of the doublet, respectively. For Na D, tb/tr=2. 
Under the curve-of-growth assumption, the optical depth r of 
a line can be expressed as: 

r(A) = Toe-(^-^°>'/(^°^/^)', (3) 

where tq and Ao are the central optical depth and central wave- 
length of the line and b is the line width. The four physical 
parameters of this model are: 

• Outflow velocity (voff), which measures the shift of line 
profile centroid relative to the galaxy's systemic veloc- 
ity (Vsys). 

• Optical depth at line center (tq) of the absorbing clouds, 
which influences the depth and intensity ratio of the two 
lines of the doublet. 

• Covering factor (C/) of the absorbing clouds, which de- 
termines the residual flux at line center 

• The Doppler width {b) of the absorption lines. 

Given the low resolution of our spectra 150 km s"' 
FWHM) and the blended nature of the profiles, there are some 
degeneracies between optical depth and covering factor which 
both influence the depth of the lines. The line width and the 
outflow velocity are more robust to degeneracies. 

In the top panel of Figure |6] we show the stacked spectra of 
Figure|5]after continuum normalization. The rest wavelengths 
of the two sodium lines are marked by vertical lines. The line 
profiles appear asymmetric, with some of the gas blue-shifted 
relative to the systemic velocity of the galaxy. We hypoth- 
esize that there is an absorption component at the systemic 
velocity that arises in the disk, and an outflow (blue-shifted) 
absorption component that arises in a galactic wind. Systemic 
and outflow absorption line components were identified in 
previous studies of individual starbursts (cf., Martin 2005), 
and in composite spectra of intermediate redshift galaxies 
dWeiner e t al. 2009) . W e co nsider the origin of the two line 
components further in 34.11 where we examine the incUna- 
tion dependence of the Na D line profile. 

We fit the Na D line profile with two absorption compo- 
nents, as shown in Figure |6] with the velocity of one compo- 
nent fixed at v'sys (red). The blue line in Figure |6] denotes the 
outflow component. We also show the single component fit in 
the top-right panel (green). Comparing the residuals of these 
two fits, we find that the two component fit is superior. To aid 
in breaking the degeneracy of the two components, we have 
also fixed the line width of the systemic component to the 
width of the He I emission line. We assume that the He I line 
profile reflects the kinematics of H II regions in the disk, and 
that the gas producing the systemic Na D absorption shares 
these kinematics. To have an idea of how the final results 
depend on this assumption, we tried two other methods: (1) 
fixing the systemic component to 1.5 times the width of the 
He I emission line; and (2) letting it vary in the range of [1- 
1.5] He I width. While method (2) gives very similar fitting 
parameters as that of fixing the width to He I, method (1) pro- 
duces larger spectral residuals and the fitted parameters show 
more scatter with galaxy physical properties. We therefore fix 
the hnewidth of the systemic component to He I. 



As noted in 32.31 there are cases where the best-fit con- 
tinuum model has stronger Na D absorption than the actual 
data. The reasons for this are still not fully understood. One 
possibility is that the young stars used in the empirical stellar 
libraries of CB08 have some excess absorption at Na D due to 
the Milky Way's ISM. However, a comparison with fully the- 
oretical models (see Appendix) suggests that such absorption, 
if present, is not very strong. Another possibility is that a com- 
bination of 40 instantaneous burst models with a range in age 
and metallicity is insufficient to represent the the true galaxy 
star formation history. Degeneracies may also limit the ability 
of our fitting code to arrive at the optimal solution. We cannot 
rule out continuum fitting errors, however, we find this expla- 
nation less likely, given the high quality of the fits elsewhere 
in the spectrum (cf., Fig.|5]l. The final possibility is that Na D 
is sometimes in emission. This is plausible because Na D is 
a resonance absorption transition. (Each absorbed photon is 
re-emitted isotropically.) Whether Na D is detected in absorp- 
tion or emission in a given spectrum depends on the relative 
geometry and velocity of the gas and stars probed by the fiber. 
We note that the emission appears slightly redshifted and it is 
seen predominantly in galaxies that are face-on and have low 
dust attenuation. This suggests that we are seeing through the 
disk to the back side of the expanding bi-polar outflow. As the 
disk becomes progressively more dusty or inclined the emis- 
sion from the far side of the disk is attenuated. Emission was 
not detected in previous studies of LIRGs and ULIRGs due 
to the dusty nature of these sources. However, Phillips ( 1993]) 
saw a clear Na D P-Cygni profile (blushifted absorption, red- 
shifted emiss ion) in the disk galax y NGC 1808. In the DEEP2 
Mg II study, IWeiner et al.l (|2009|) found Mg II emission in a 
subset of the galaxies. The emission tends to be in the bluest 
of the galaxies. In some of them it is visible in the individual 
spectra, not just the stack. They don't have inclination infor- 
mation for the high-z galaxies, but it is natural to think that 
bluer galaxies are lower inclination on average. 

We include the objects with emission Na D residuals in the 
stacking and accommodate the possibility of Na D emission 
by allowing the systematic component in our line profile fit 
to be in either absorption or emission. We note that in some 
cases a substantial systemic emission component is fit, even 
when the line profile appears to be completely in absorption 
(cf., the first panel of Figure[8]l. In such cases, large residuals 
result from imposing the requirement that the systemic com- 
ponent be in absorption (note the strong asymmetry of the line 
profile). We believe these fits are robust and represent real 
cases of superimposed Na D emission from the disk and ab- 
sorption from the outflow. Figure|9]shows very smooth trends 
with galaxy physical properties as the systemic component 
transitions from absorption to emission, suggesting that our 
code generally does a good job of distinguishing the outflow 
and systemic components whether the latter is in absorption 
or emission. We have therefore chosen to allow the systemic 
component to be in emission; for simplicity, we implement 
this by allowing for negative covering factors. 

4. PROBING GALACTIC WINDS USING STACKED 
SPECTRA 

The dependence of the outflow properties on galaxy physi- 
cal properties is of great interest because it can help clarify the 
influence of SFR and stellar mass on the wind kinematics. We 
can also learn which galaxies host outflows and thereby gain 
insight into the origin of metals in the IGM. These empirical 
results can, in turn, be incorporated into numerical Simula- 
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Fig. 7. — Continuum-normalized spectra binned by galaxy inclination. Dif- 
ferent inclination bins are indicated by different colors, with the galaxies clos- 
est to face-on in red and the galaxies closest to edge-on in black. 



tions of galaxy formation and evolution. 

The data analysis method described in §3 is applied into 
our star forming galaxy sample in this section. We explore 
the effect of changes in galaxy inclination in ^4.11 and ex- 
amine how the outflow and systemic components of the ISM 
Na D absorption evolve with various galaxy physical proper- 
ties in ^4.21 This is the first time outflow properties have been 
studied with large complete sample in the local universe that 
includes normal star forming galaxies as well as starbursts. 

4. 1 . Trends with Galaxy Inclination 

It has been suggested that the outflow component of the 
Na D absorption arises from ambient interstellar material that 
has been entrained and accelerated along the minor axi s of the 
galax y by a hot starburst-driven superwind (Hec kman et alj 
120001) . while the systemic component of the absorption arises 
in the disk. Based on this simple picture, we expect the two 
Na D components to have strong - and opposite - variations 
with galaxy inclination. We test our assumptions about the 
absorber geometry using spectra stacked in bins of galaxy in- 
clination. 

Galactic winds are ubiquitous in galaxies with Esfr > 
O.IMq yr"' kpc~'. To test our picture of the wind outflow 
geometry, we therefore apply this Ssfr cut to our 'Sample 
A' galaxies, to select sources likely to host outflows. In Fig- 
ure|2]we compare the interstellar Na D absorption profiles for 
these galaxies in 8 different inclination bins. Three features 
are immediately obvious. First, the line profile depth at the 
rest wavelength of the doublet increases dramatically as incli- 
nation increases and galaxies are viewed more nearly edge-on. 
Second, the line center shifts to bluer wavelengths as inclina- 
tion decreases and galaxies are viewed more nearly face-on. 
Third, the lines are not saturated — in which case the dou- 
blet ratio would be closer to 1:1 - nor are the profiles black 
at the line center, which would indicate complete covering of 
the source by the absorber 

Figure [8] illustrates our two-component absorption line fits 
to composite spectra in different inclination bins. The red 
lines show the component fixed at Vjys, the blue lines show 
the outflow component, and the green lines show the com- 
bined fit. The systemic component (red) evolves markedly 
from emission at low inclinations to become the dominant ab- 



sorption component at high inclinations. The outflow compo- 
nent (blue) evolves strongly in the opposite sense, dominating 
the line profile in face-on systems, but becoming a minor con- 
tributor to the absorption in edge-on galaxies. 

The equivalent width (EW), velocity (voff), line width (b), 
covering factor (Cf) and optical depth at line center (tq) of 
the two absorption line components are shown as a function 
of inclination angle in Figure |9l At / < 60°, the outflow ve- 
locity depends only weakly on inclination, and has a median 
value of ^-^140 km s"'. Above / = 60°, the velocity drops 
precipitously. The line width of the outflow component de- 
creases by 40% as galaxies go from face-on to edge-on, but 
shows considerable scatter. The EW of the outflow compo- 
nent changes smoothly and nearly linearly from EW=0.7 A 
to at / ^ 0° to nearly 0.1 at / = 90°. This trend appears to 
be driven primarily by a decrease in the gas covering factor 
The optical depth of the outflow component increases sharply 
above / = 60°. This may indicate that the wide-angle out- 
flow of dusty material has a more transparent and fast-moving 
"core" (/ < 60°) and a surrounding sheath of optically-thick 
slower-moving entrained material (/ > 60°) - sort of a "turbu- 
lent boundary layer". However, the outflow component makes 
a relatively small contribution to the total EW in this regime, 
so some caution in interpreting the result is advised. 

The systemic absorption component changes with inclina- 
tion in the opposite sense of the outflow component: the line 
width increases with inclination by ^ 20% and the EW in- 
creases, from emission at / ^ 0° to 0.7 A of absorption at 
/ ~ 90°. The lack of zero-velocity ISM absorpti on component 
at low inclinations is consistent with the result of Weiner et alJ 
(|2009), in which they found that in the bluest subset, which 
are likely to be galaxies viewed face-on, there was no com- 
ponent of absorption at systemic velocity. Again, most of 
this change can be attributed to the covering factor; the opti- 
cal depth shows little trend with inclination, and considerable 
scatter The scatter in Figure |9] represents a lower bound to 
the uncertainities in the fits. 

The trends shown in Figure|9]are in agreement with our ba- 
sic hypothesis that Na D absorption arises in both the galaxy 
disk and in an outflow directed along the galaxy minor axis. 
Notably, the strength of the systemic component increases 
as galaxies are viewed more nearly edge-on and a greater 
path length through the disk is probed. This leads us to con- 
clude that the systemic component arises in the warm neutral 
medium of the disk. These disk clouds could either be in the 
central star forming region, in which case they must be too 
massive or too pressure-confined to be affected by feedback 
from young star clusters, or the clouds could be in the outer 
more quiescent part of the disk which is probed by the fiber at 
high galaxy inclinations. 

The outflow absorption component covers a wide angle, but 
seems stronger within ^ 60 ° of the disk rotation axis; it likely 
arises from a large-scale galactic wind, like that seen in Fig- 
ure [3] Notably, the velocity is roughly constant at / < 60°, 
as expected if the outflow is radially directed. At / = 60°, 
the EW of the disk-like and outflow components are roughly 
equal. This is consistent with our understanding of the out- 
flow and disk geometry: at / < 60°, the ISM Na D absorp- 
tion is dominated by outflow, while the disk-like component 
becomes important at / > 60°. Based on these results, we 
will split "Sample A" into two sub-samples, "Sample B", with 
/ > 60°, and "Sample C" with / < 60°. There are 60,21 1 and 
80,414 galaxies in each sample, respectively. 
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Fig. 8. — Absorption line profile fits to composite spectra in different inclination bins. The black lines show the continuum normalized composite spectra; the 
red line shows the absorption line component fixed to the galaxy systemic velocity; the blue fine shows the outflow component; and the green line shows the 
combined fit. Note the appearance of Na D emission at the systemic velocity in the low inclination bins. 



4.2. Trends with Galaxy Physical Properties 

In this section we examine how the properties of the two 
Na D absorption components (as derived from our constrained 
two-component fit) vary as a function of galaxy physical prop- 
erties. "Sample B" and "Sample C" are used to study the disk- 
like and outflow components, respectively. Our reasoning is 
that the absorbers in the disk are best studied in galaxies that 
are more nearly edge-on, while absorbers in the outflow are 
best studied in galaxies that are more nearly face-on. 

We consider each galaxy physical parameter independently 
and use the adaptive binning approach described in ^3.11 to 
construct stacked spectra. We do not consider the full param- 
eter range, but rather we define a lower boundary based where 
the fraction of strong Na D absorbers (EQ > 0.8 A) begins 
to rise steeply (see the blue lines in Fig.|4]i. The continuum- 
normalized spectra for eight different Ssfr bins from "Sample 
C" are shown in Figure [lO] as an example. The strength of the 
absorption increases dramatically with increasing Ssfr- Fig- 
urefTllshows the two-component fits of these eight spectra. 

We plot the parameters derived from our absorption line fits 
(voff, b, EW, Cf, To) as a function of galaxy physical properties 
(SsFR, Ay, SSFR, M*) for Samples B and C in FiguresfT2land 
[T3] We use these plots to help discern the galaxy physical pa- 
rameters that have the greatest influence on the warm neutral 
medium in and outside of galaxy disks. However, the strong 



correlations among the galaxy physical parameters illustrated 
in Figure [H make this exercise challenging. We also consid- 
ered correlations with SFR and S^, (both measured in the 
fiber aperature). Trends with these parameters are very sim- 
ilar to those found with Ssfr and M*, but typically weaker 
and/or nosier For clarity we have not included SFR and Y.m, 
in our subsequent analysis. 

For the systematic component (Fig. ITS]) , we note the fol- 
lowing: 

• Na D EW increases strongly and nearly Unearly with 
log(I]sFR), log(M*), and Ay from negative values 
(emission) to ^ 1 .0 A. We note that the line EW is sensi- 
tive to the velocity spread of the gas, the cloud covering 
factor and the optical depth. 

• The line width (which is tied to the width of the He I 
line) also increases smoothly with all of the stacking 
parameters except SSFR. The strongest trend seen is a 
60% increase in b with stellar mass. Since Sample B 
galaxies are highly inclined, the line width probably re- 
flects the rotation speed of the inner disk. The apparent 
correlation between b and Ssfr, Ay most likely results 
from the correlation of M^, with these parameters, as 
show in Figure [T] 

• The covering factor changes steeply, from negative val- 
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Fig. 9. — Derived parameters of the outflow and systemic Na D absoiption components as a function of inclination angle. The blue diamonds represent the 
outflow component and the red circles represent the systemic component. The parameters shown include the velocity of the outflow component (I'ofi), the Une 
width (b), the covering factor (Cf), the optical depth at line center (tq) and the Hne equivalent width (EW). 
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Fig. 10. — Continuum normalized composite spectra of galaxies in 8 dif- 
ferent SspR bins drawn from Sample C (i < 60°). As SsFR increases, the 
absorption becomes stronger SgpR increases from red to purple. 



line profiles. 

• The optical depth ranges from T()=0 to ^ 1.5, and it does 
not depend strongly on any of the stacking parameters. 
In general, the optical depth is the least well constrained 
parameter that we measure. 

Examining the changes in b, Cf, and tq, we hypothesize 
that the changes in EW are driven primarily by the covering 
factor with an small additional contribution due changes in 
the line width. It is difficult to determine whether Esfr, Ay, 
or are the primary driver of the trends. SSFR is clearly not 
an important parameter. We address the iss ue of correlations 
between the various physical parameters in 34.31 

In Figure [T3] we explore the dependence of outflow Na D 
component on various galaxy properties. We note the follow- 
ing: 

• The Na D EW increases with Ssfr, M^, and Ay (similar 
to the disk component). 

• The line width appears to scale strongly with Esfr, 
and Ay. The tightest correlation appears to be with : 



ues (Na D emission) to C/ = 0. 15. The very small Na D 
covering factor may be related to the low filling factor 
of neutral gas in galactic disks. The trends are similar 
to those identified with the EW, but noisier due degen- 
eracies with To exacerbated by the blended nature of our 



• The covering factor is very low (Cf ^ 0.04-0.12) and 
changes by a factor of three, at most. The strongest 
correlation is with Ay and M* . 
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Fig. 1 1. — Absorption line profile fits to composite spectra in different Ssfr bins. (The spectra are the same as those shown in Figure [To]) The black lines 
show the continuum-normalized spectrum; the red line is the systemic (disk) absorption component; the blue line is the outflow component; and the green lines 
is the combined model. The range in log SspR is given in the lower right of each panel. 



• The lines range from unsaturated to mildly saturated 
(ro ^ 3) and tq displays a pronounced inverse corre- 
lation with Ay, and M,. 

• The outflow velocity does not show any strong trend 
with galaxy physical properties. A very weak and noisy 
correlation is evident with Ssfr and possibly with M^, . 

• No properties of the outflow component are well corre- 
lated with the specific star formation rate. 

One possible interpretation of the observed trends in C/ and 
To is that in galaxies with low dust attenuation Na D only sur- 
vives in very dense clouds in the outflow. These clouds are 
optically thick and comparatively rare, so the average tq is 
high and the Na D covering factor is very low. In galaxies 
with higher dust-to-gas ratios, Na D is shielded from ioniz- 
ing radiation in clouds with lower optical depths. Since these 
clouds are more numerous, the Na D covering factor is higher 
and the average optical depth is lower Thus, the mean optical 
depth decreases with Ay while the covering factor increases, 
as observed. 

The lack of trends with SSFR is somewhat surprising since 
it might be expected that outflows depend on both the SFR and 
the potential well depth. One explanation is the strange shape 



of the SSFR vs. Ay correlation (FigdJ. Galaxies with both 
low and high SSFR have low dust attenuation. This arises 
because low SSFR galaxies are predominately gas-poor mas- 
sive galaxies with high dust-to-gas ratios, while high SSFR 
galaxies are predominantly gas-rich dwarfs with low metallic- 
ities and dust-to-gas ratios. At intermediate values of SSFR, 
galaxies display a wide range in dust attenuation. Clearly, 
dust shielding plays a key role in the detectablity of outflows 
using Na 1. 

The EW of the outflow component will depend on b, Cf, 
and To, although the influence of tq is minimal for galaxies 
with To > 1. The EW depends linearly on b and Cf, but the 
dynamic range probed by our sample is greater in C/ (factor of 
^ 6 vs. factor of ~ 3.) The EW thus primarily reflect changes 
in the covering factor. Since our measurement of the EW is 
less noisy than our measurement of Cf, we use the EW in ex- 
amining tends with galaxy physical properties further below. 

4.3. Isolating the Primary Drivers of the Observed 
Correlations 

In figures [12] and [T3] we identified Ay, Ssfr, and as 
important physical parameters controlling the Na D line pro- 
file. However, due to strong correlations among these param- 
eters (Fig. [U, it is difficult to isolate the primary drivers of 
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the trends we observe. The correlations are sufficiently strong 
that it is not possible, for example, to look for trends with 
EsFR at fixed M*. A bin with a small range in stellar mass 
would also have a small range in Ssfr- To isolate the most 
important physical parameters, we adopt the following ap- 
proach. We take each bin in M» (for example) and further di- 
vide it into three equal sub-bins, sorting the galaxies by Esfr- 
We are then able to examine trends between Na D absorption 
and in low, medium, and high Ssfr bins (where the exact 
division between low, medium, and high changes with M,). 
If is the primary driver of the trend, then no difference 
is expected between the sub-bins. Conversely, if Ssfr is the 
primary driver of the trend, then the three sub-bins will be 
strongly offset from one another at each value of M» . 

Figure [14] shows the Na D EW of the systemic component 
as a function of M», Ssfr, and Ay for galaxies in Sample B 
with log M^JMq > 10.3. In each panel the sample is split 
into three sub-bins according to the galaxy parameter labeled 
in the top-left corner. The black, red, and blue points indicate 
the low, medium, and high sub-bins respectively. The dom- 
inant driver of the observed trend (the parameter displaying 



the smallest offsets between the sub-bins) appears to be Ay. 
This suggests that high dust attenuation is important to keep 
Sodium in the disk neutral. A small residual correlation is 
evident with Ssfr and M*. A correlation between Ssfr and 
the Na D covering factor is expected if Esfr scales with the 
filling factor of molecular clouds, as suggested in Bigie l et akl 
(2008). The residual trend of EW with stellar mass is due 
to the larger velocity spread of the line in disk galaxies with 
larger rotation speeds. 

Figure [15] shows the Na D EW of the outflow component 
as a function of M,, Ssfr, and Ay for galaxies in Sample C 
with log(M,/M0) > 10.3. As in Figure [T4l the galaxies are 
split into three sub-bins according to various physical proper- 
ties. The plots with stellar mass on the x-axis show very large 
discrepancies between high and low sub-bins in Ssfr and Ay, 
therefore, we infer that is the least important physical pa- 
rameter in determining the EW. The most important parameter 
influencing the line EW appears to be Ssfr followed closely 
by Ay. The physical explanation for these correlations is 
straightforward: star formation surface density determines the 
amount of material lofted above the disk and the dust attenu- 
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ation influences the survival of Na I in the clouds. 

A similar analysis applied to the outflow velocity (v'off) in- 
stead of the EW does not yield any insight into the domi- 
nant physical parameter determining the outflow speed. We 
note that the dynamic range in Voff is very small (120 - 160 
km s"'), and this makes identification of trends particularly 
challenging. It is also important to keep in mind that the range 
in galaxy physical properties probed by our study is quite 
small - typically of order 1 dex. For instance, we are look- 
ing for trends over a range in stellar mass of log(M»/M0) = 
10.3 -11.3, whereas the full galaxy population spans over six 
orders of magnitude in stellar mass. This limitation is im- 
posed by our need to study objects where Na D absorption is 
common (Fig.|4|i. 

In spite of these limitations, we can use our results to in- 
vestigate a basic phy sical model of the outflow. Following 
iMartin & Bouchel (|2009|) we suggest that the velocity of the 
line centroid reflects the speed of the swept-up shell of inter- 
stellar gas at the point whe re it blows out of the dis k. Combin- 
ing equations 3 and 6 from lStrickland etaTI (120041) . for the ve- 
locity of a sheU of gas at blow-out, we find Vos oc (Ssfr/po)'^"^ 



where po is the gas density. Assuming po oc Y^gas/H^, where 

H, is the scale height of the disk, and E^^, cx Sgp^"^ (iKennicutj 

fT998bl) . we find V'off oc E^pj^//; This is consistent with the 
very shallow scaling we see between outflow velocity and 

SsFR- 

Once the cool entrained gas has escaped the disk, it may 
be accelerated further by the ram pressure of the hot wind, 
or by radiation pressure on dust grains. T his halo gas likely 
contributes to the wings of the line profile (iMartin & Bouchel 
l2009h . In Figure [m we examine the physical drivers of the 
line width of the Na D outflow component. We see that the 
line width is sensitive to all three parameters: M*,S]sfr, and 
Ay. The latter two are intuitive: cloud acceleration is likely 
related to Esfr, while dust attenuation is important to keep the 
gas neutral. The positive correlation between line width and 
stellar mass is more puzzling. It may be that clouds with ve- 
locities above the halo escape velocity quickly dissipate after 
leaving the pressure-confining hot medium of the halo. The 
larger line widths in more massive galaxies may therefore re- 
flect the larger size of the hot halo. 
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Fig. 14. — The dependence of the Na D EW of the systemic component 
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4.4. Comparison with Previous Work 

A correlation between outflow velocity and galaxy SFR 
and rotation s peed was identified bv lRupke et al.l (l2005b t) and 
iMartinI (l2005h in starburst galaxies. In FigurefTTlwe compare 
our data (cyan crosses) to these results. For consistency with 
the literature data, we use global (i.e., aperture-corrected) 
SFRs. Green asterisks show 1 Jy ULlRGs from Rupke et al] 
( 120021) : red squares indicat e ULIRGs from Martin (2005^ 
blue circles are LIRGs from lHeckman et aTl (2000); pink di- 
amonds are z 1.4 galaxies from ("Weiner et al. 2009); and 
black triangles are dwarf starbursts from Schwartz & Martin 
(12004'). All the data are corrected to a Kroupa IMF. 

The velocities of our composite spectra show good 
agreement with the velocities der ived from ind i vidual 
galaxies in the same S FR range dHeckman et all llOOOl; 
[Schwartz & Martini |2004l) . However, fli e trend wifli S FR is 
weaker than the correlation suggested by 'Martin' (2005, black 
line). We note that slope of this correlation is strongly influ- 
enced by three dwarf galaxies with very low SFRs. The large 
scatter on Voff at fixed SFR for the individual measurements 
may be due both to intrinsic variations and inclination effects. 

Based on a sample of high stellar mass galaxies between 
z ~ 1.4 and z ~ 1, [Rubin et al.. (.2009.) suggest that outflow 



Fig. 15. — The dependence of the Na D EW of the outflow component 
on galaxy physical properties. As in Figure [14] the sample in each panel is 
divided into three sub-bins by the parameter marked in the top-left comer 
(black=low, red=medium, blue=high). 



absorption strength (measured from Mg II AA2796,2803) is 
more closely associated with SFR than with M,, and it does 
not increase with increasing SSFR. Although they have not 
detected a strong dependence of outflow absorption strength 
on EsFR, they have by no means ruled it out due to the spec- 
ti-um S/N limit. 

5. SUMMARY 

In this paper, we study intersteflar Na I "D" AA5890,5896 
absorption in a large sample of star forming galaxies drawn 
from SDSS DR7. At the low SFRs probed by our sample, 
cool stars with prominent stellar Na D lines make a signifi- 
cant contribution to the integrated light. We account for this 
by modeling the full galaxy continuum using the CB08 stel- 
lar population synthesis models. We find a high incidence of 
strong ISM Na D absorption lines (EW > 0.8 A) in galaxies 
that are massive, heavily dust attenuated, and that have high 
SsFR(Fig.|4]) 

We use the spectral stacking technique to increase the 
signal-to-noise ratio of the spectra and the completeness of 
the sample and stack spectra in various bins of galaxy physi- 
cal properties. In the continuum-normalized spectra we iden- 
tify two interstellar Na D absorption components — one at 
the systemic velocity, which is most pronounced in edge-on 
galaxies, and one outflow (blue-shifted) component which is 
dominant in face-on galaxies. We use two-component absorp- 
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Fig. 16. — The dependence of the line width of the outflow component 
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tion line profile fits to measure the outflow velocity, the ab- 
sorption line EW, the line width, covering factor, and the op- 
tical depth. As highlighted in the appendix, our results are 
robust with respect to the modeling of the stellar continuum. 
We find that: 

• The ISM Na D absorption arises from cool gas in the 
disk, and from cool gas entrained in a galactic wind that 
is perpendicular to the disk and has an opening angle of 
-60°. 

• For the systemic (disk) component, the Na D EW de- 
pends most strongly on the dust attenuation. Dust at- 
tenuation is undoubtedly important in shielding Na 1 
from the ionizing radiation of OB stars. At fixed Ay, 
there is an additional dependence on Ssfr- Galaxies 
with higher Ssfr probably have a higher filling factor 
of cold clouds and thus a higher Na D covering fraction. 

• For the outflow component the Na D EW depends 
strongly on Esfr and secondarily on Ay. We hypothe- 
size that the star formation surface density determines 
the amount of material lofted above the disk, and the 
dust attenuation influences the survival of Na 1 in the 
clouds. 

• The covering factor of Na I in the outflow increases 
with increasing dust attenuation while the optical depth 
of the absorption decreases. This suggests that Na D is 
able to survive in clouds with lower column density as 
the dust-to-gas ratio increases. 

• The outflow velocity (line centroid) does not depend 
strongly on any galaxy physical parameters over the 
limited dynamic range in physical properties probed by 
our study. There is some evidence for a very shallow 
trend, Voff oc Sgpj^, which is consistent with theoretical 
expectations for the velocity of a swept-up shell of gas 
at the point where it blows out of the disk. 

• The line width of the outflow component is sensitive 
to EsFR, Ay, and M,. These parameters are likely to 
influence the acceleration of halo clouds (by ram or ra- 
diation pressure), the ability of clouds to remain neu- 
tral, and the length of time a cloud can remain pressure 
confined. 

The question of how galactic scale outflows are influenced 
by AGNs is also of great interest. In a companion paper, we 
will study the Na D absorption properties in SDSS galaxies 
hosting type 2 AGNs using methods similar to those described 
here. 

The Na D ISM absorption line provides a convenient probe 
of galactic winds due to its location in the optical. However, 
the low EW width of the line in typical galaxy spectra and 
the strength of the Na D stellar feature make interstellar Na D 
absorption a challenge to measure accurately. The strong sen- 
sitivity of Na 1 to ionization by OB stars is also a very se- 
rious concern, in particular for studies that aim to compute 
mass outflow rates from Na 1 column densities and veloci- 
ties. Higher redshift studies that make use of higher ioniza- 
tion transitions such as Mg 11 A2796,2803 may ultimately be 
able to provide more robust constraints over a larger range in 
galaxy physical properties. 

Theoretical models of galaxy evolution have begun to in- 
corporate galactic winds, using various ad-hoc prescriptions 
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based on our knowledge of the cool gas. However, the ma- 
jority the energy and newly-synthesized metals in the outflow 
resides in the hot X-ray emitting phase of the wind (D. K. 
Strickland & D. C. Dinge, 2010, in prep.). To truly under- 
stand the impact of galactic winds on the evolution of galax- 
ies and the intergalactic medium, it is crucial to measure the 
chemical composition and velocity of the hot gas. This will 
require a high sensitivity X-ray imaging spectrometer such as 
the planned International X-ray Observatory. 
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Fig. 18. — This figure is the same as Figure[5] but we add the GD05 stellar continuum model (blue line) for comparison. The CB08 model provides a better fit 
to the stellar population. 
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APPENDIX 

Here we compare the stellar continuum of lGonzalez Delgado et alj (120051 hereafter GD05) and CB08 stellar population models 
to check whether our results are model dependent. 

In §2.1, we noted that stellar Na D appears to be over-estimated in our continuum fits, resulting in emission at the systemic 
velocity in the continuum-normalized spectra. The libraries of stellar spectra used by CB08 are empirical, so Na D could be 
contaminated by Milky Way absorption. One important question is to what extent our results are affected by the stellar population 
models used to estimate the continuum. In order to answer this question, we repeat all our analysis using the models of GD05 
which use fully theoretical stellar libraries. Figure [18] shows both CB08 and GB05 model fits to the stacked spectrum shown in 
Figure |5] The CB08 model is shown in red and the GD05 model in blue. The CB08 model clearly provides a better fit of the the 
full spectrum and the region around the Mg 1 triplet. 

Using the GD05 theoretical models, all the trends shown in Figs. |9] [12] and [13] remain. On average, the outflow velocity 
decreases by ~ 40-60 km s~' ; the Na D EW increases by ^ 0.3 A, for both disk-like and outflow components; and the covering 
factors for the two components increase by 10%. The line widths for the systematic component are the same as that from CB08, 
since they are bound to He 1. The width of outflow component increases by ^ 25 - 30 km s"' and more scatter is introduced in the 
trends with galaxy physical properties. The largest change caused by the using different stellar population models is the optical 
depth tq. However, this parameter is not very strongly constrained in either case. In summary, the CB08 models provide a better 
overall fit to the stellar continuum. However, our basic results are not sensitive to the stellar population models used. 



